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In the Nastrom-Gage spectrum of atmospheric turbulence we observe a k 3 energy spec- 
trum that transitions into a A: -5 / 3 spectrum, with increasing wavenumber k. The tran- 
sition occurs near a transition wavenumber k tl located near the Rossby deformation 
wavenumber kn- The Tung-Orlando theory interprets this spectrum as a double down- 
\mJ scale cascade of potential enstrophy and energy, from large scales to small scales, in which 

the downscale potential enstrophy cascade coexists with the downscale energy cascade 
over the same length-scale range. We show that, in a temperature forced two-layer quasi- 
geostrophic model, the rates with which potential enstrophy and energy are injected 
place the transition wavenumber k t near k^. We also show that if the potential energy 
■^j- dominates the kinetic energy in the forcing range, then the Ekman term suppresses the 

upscale cascading potential enstrophy more than it suppresses the upscale cascading en- 
ergy, a behavior contrary to what occurs in two-dimensional turbulence. As a result, the 
ratio rj/e of injected potential enstrophy over injected energy, in the downscale direction, 
decreases, thereby tending to decrease the transition wavenumber k t further. Using a 
random Gaussian forcing model, we reach the same conclusion, under the modeling as- 
sumption that the asymmetric Ekman term predominantly suppresses the bottom layer 
forcing, thereby disregarding a possible entanglement between the Ekman term and the 
t-h nonlinear interlayer interaction. Based on these results, we argue that the Tung-Orlando 

theory can account for the approximate coincidence between k t and k^. We also identify 
certain open questions that require further investigation via numerical simulations. 

% 



1. Introduction 

Quasi-geostrophic models capture the dynamics of the atmosphere at planetary scales 
greater than 100km, in order of magnitude. They are based on the assumptions of rapid 
rotation and small vertical thickness, both of which become pronounced features of the 
dynamics of atmospheric motion at increasingly larger length scales. The simplest quasi- 
geostrophic model is the two-layer model, in which we have two layers of two-dimensional 
vorticity-streamfunction equations, coupled by a temperature equation, situated in a mid- 
layer between the vorticity layers. Obviously, the two-layer model can be generalized by 
adding more layers of vorticity-streamfunction equations interlaced with temperature 
equation mid-layers. In the limit of an infinite number of layers, we converge to the full 
quasi-geostrophic model. 
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Until recently, it was assumed that quasi-geostrophic turbulence has the same dynam- 
ical behaviour as two-dimensional turbulence, where, according to the theory of |KraichT 
^p6^,|L^ithlp68l), 



nan 



and Batchelor (1969), there is a downscale enstrophy cascade 
and an upscale inverse energy cascade. This assumption follows from |Charney| ( 1971 1 
who claimed that there is an "isomorphism" between quasi-geostrophic turbulence and 
two-dimensional turbulence. More recently, key differences between these two models 



were noted by Tung & Welch (20011, Tung & Orlando (20036), and Gkioulekas & Tung 



(2007a). The most remarkable difference was highlighted in the numerical simulation of 
the two-layer quasi-geostrophic model by Tung & Orlando (2003a), which produced an 
energy spectrum that scales as fc~ 3 initially, and with increasing wavenumber k, tran- 
sitions to A; -5 / 3 scaling. This is consistent with the observed energy spectrum of the 



atmosphere, as was first measured by Nastrom & Gage (1984) and Gage & Nastrom 



(1986), but it is inconsistent with our conventional understanding of the dynamical be- 
havior of two-dimensional turbulence, as described by the theory of Kraichnan (1967), 



Leith (1968), and Batchelor (1969). Tung fc Orlando (2003a I showed that their sim- 



ulation produced a downscale enstrophy cascade that co-existed with a downscale en- 
ergy cascade, with both energy and enstrophy injected by baroclinic instability at small 
wavenumbers, and dissipated at large wavenumbers. If 77 is the enstrophy flux and e is 
the energy flux as s ociated with these coexisting cascades, then, by dimensional analysis, 



Tung & Orlando (2003 a 



argue that the transition from k~ A scaling to fc~ 5 / 3 scaling 
should occur at a transition wavenumber k t ~ (77/e) 1 / 2 , and it does. 

This result bolstered the Tung-Orlando theory that interpreted the Nastrom-Gage 
spectrum as a coexisting downscale potential enstrophy cascade and downscale energy 
cascade, both spanning a comparable range of length scales. It should be noted that it 
was not the intention of Tung & Orlando (2003a) to claim that the entire Nastrom-Gage 
spectrum can be explained via the two-layer quasi-geostrophic model. The point of the 
simulation was to demonstrate that it is possible to have coexisting downscale potential 
enstrophy and energy cascades, even in models as close to two-dimensional turbulence 
as the two-layer quasi-geostrophic model. This possibility is bound to become even more 
favorable under models that are further away from the two-dimensional approximation, 
such as the multi-layer quasi-geostrophic model or the three-dimensional stratified turbu- 
lence model. Gravity waves can also play a helpful role in facilitating coexisting cascades, 
as discussed further in section 5. 

In a subsequent paper, Smith (2004) criticized Tung & Orlando (2003a) on the follow- 
ing grounds: First, using his "HVK scale" estimates, he claimed that the Tung-Orlando 
numerical simulation is not well-resolved and that therefore the A: -5 / 3 part of the Tung- 
Orlando energy spectrum is a bottleneck instead of being indicative of a real cascade. 
In connection with this clai m, |Smith ( 2004 1 criticized the use of a resolution-dependent 
hyperdiffusion coefficient by Tung & Orlando (2003a). Second, that in two-dimensional 



turbulence it is not possible for the downscale energy flux to be large enough to cre- 
ate a gap between the transition wavenumber k t and the dissipation wavenumber fe<j. In 
response, Tung (2004) noted that: (a) Since, the diagnostics in the Tung-Orlando sim- 
ulation indicate a downscale enstrophy flux 77 and a downscale energy flux e that have 
magnitudes consistent with the location of the transition wavenumber kt in the simu- 
lation's energy spectrum, it is very unlikely that the transition is caused by an energy 
bottleneck, as argued by Smith (2004 1, (b) The use of a resolution-dependent hyper- 
viscosity coefficient is intended to model the anomalous energy dissipation sink at small 
scales, originating from three-dimensional dynamics, by controlling the downscale energy 
dissipation rate, (c) Smith's "HVK scale" argument, which was used to argue that the 
Tung-Orlando simulation is not well-resolved, has various flaws, discussed in detail by 
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Tung (2004), that render it inconclusive. Nevertheless, Smith's claim, that kt will coin- 
cide with the dissipation scale A 1 ,/, can be shown to hold, for the case of two-dimensional 
turbulence, via a corrected proof given by Gkioulekas k Tung (20056). However, as 
was shown by |Gkioulekas k Tung (2007a), this result is not necessarily generalizable 
to quasi-geostrophic models. Thus, Smith (2004) did not establish the claim that "an 
inertial range transition is not possible in quasi-geostrophic models" , and the theoretical 
problem remains open, (d) Since the diagnostics of the Tung-Orlando simulation indicate 
that the downscale energy dissipation rate balances the rate with which energy is sent 
downscalc from the forcing range, the simulation is sufficiently well-resolved to prevent a 
bottleneck-type energy pile up at small scales, over a time-scale longer than the runtime 
of the simulation. 

Be that as it may, the underlying theoretical question, implied by Smith (2004), re- 
mained open: How can the downscale energy flux e be large enough to yield a gap 
between kt and kd, when that is not possible in two-dimensional turbulence? An even 
deeper question also demanded further understanding: how is it possible for two down- 
scale cascades to coexist? These lingering questions generated skepticism towards the 
Tung-Orlando theory, which is why we were prompted to investigate them at greater 
depth. In |Gkioulekas k Tung| ( [2005a|6[ ), we augmented the Tung-Orlando theory by not- 
ing that even in two-dimensional turbulence there is a small amount of energy cascading 
from small to large wavenumbers, as long as the viscosity coefficient of the small-scale 
dissipation term is non-zero. We have proposed that this small "energy leak" should 
be viewed as a downscale energy cascade that coexists with the dominant downscale 
enstrophy cascade. To support this theory, in |Gkioulekas k "T ung (200 5a|6 ) we noted 
that the triad interactions responsible for the enstrophy cascade are independent from 
those responsible for the downscale energy cascade. This is, in fact, an immediate but 
unstated consequence of the original argument by Kraichnan (19671, as noted in section 



3.2 of Gkioulekas k Tung (2005a). We have also theorized that the two cascades can 
be viewed as two independent homogeneous solutions of the governing statistical theory 
that can be linearly superposed on each other. |Davidson| ( |2008[ ) confirmed the validity of 
the linear superposition principle for the 3rd-order structure functions, thereby adding 
further detail to a corresponding proof sketch given in section 3.1 of Gkioulekas & Tung 



(2005a). Some of the details of my statistical theory of two-dimensional turbulence was 



given in Gkioulekas (2008) and Gkioulekas (2010), and further development of this theory 



is currently in progress. 

To elaborate further, our claim is that the energy spectrum of the downscale cascade 
is given by the linear combination of a dominant k~ 3 term, arising from the dominant 
downscale enstrophy cascade, and a subdominant fc~ 5 / 3 term, arising from the hidden 
downscale energy cascade, which allows, in principle, a transition from the —3 slope to 
the —5/3 slope. In linearly dissipated two-dimensional turbulence, this transition is not 
expected to be realized, because an upper bound on the energy flux forces the transition 
scale kt to be greater than the dissipation scale kd of the enstrophy cascade. If HE{k) is 
the energy flux from the (0, k) wavenumber interval to the (k, +oo) interval and Hoik) 
is the enstrophy flux from (0, k) to (fc, +oo), then this flux inequality reads k 2 IlE(k) — 
^c{k) < 0, for all wavenumbers k not in the forcing range. The origin of this inequality 
is the relationship D E {k) = k~ 2 D G (k) between the energy dissipation spectrum D E {k) 
and the enstrophy dissipation spectrum DQ(k). Thus, with increasing wavenumber k, 
the proportion of the energy dissipation rate relative to the enstrophy dissipation rate 
vanishes rapidly with k~ 2 , and this rapid decrease is the reason why the downscale energy 
cascade coexisting with the downscale enstrophy cascade cannot be seen in the energy 
spectrum. However, the subtle point that deserves to be stressed here is that the proof of 
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the inequality involves both the linearity of the dissipation terms and the twin conservation 
laws (of enstrophy and energy), and is not inherent solely to the twin conservation laws. 
With a nonlinear dissipation term, that could result from neglected three-dimensional 
effects, the flux inequality could be violated within the inertial range. Since the transition 
wavenumber fc t is expected to be approximately equal to the wavenumber where the flux 
inequality becomes an equation, an inertial range violation of the flux inequality would 
give kt < k d . 

in the quasi-gcostrophic two-layer model, 



As was shown by Gkioulekas & Tung ( 2007a 



the relationship between DE{k) and D G {k) becomes quite complicated, so it may be pos- 
sible to violate the flux inequality, thus resulting in a significant separation between kt 
and kd- If that occurs, we can expect fc -5 / 3 scaling in the gap created between kt and 
kd- In Gkioulekas & Tung (2007a), we have identified asymmetric dissipation as the only 



mechanism that can break the flux inequality. By asymmetric dissipation we mean that 
the dissipation operators acting on the top and bottom velocity equations are different: in- 
deed, in the two-layer quasi-geostrophic model there is an Ekman dissipation term acting 
at large scales at the bottom layer but not at the top layer. Since the small-scale hyper- 
diffusion is not physically inherent in the quasi-geostrophic dynamics, there is no physical 
reason to prefer symmetric over asymmetric hyperdiffusion either. Unfortunately, there 
are still many open questions concerning the theory of the flux inequality. Consequently, 
the numerical results of Tung & Orlando (2003a) notwithstanding, there is still some 
uncertainty on whether the two-layer model can have a robust energy dissipation sink 
that can break the flux inequality in the inertial range. 

It should be noted that when the same dissipation operator is used on both layers, it 
can be proved that the flux inequality is satisfied for all wavenumbers not in the forcing 
range. For that case, the prediction of Charney ( 1971 ), that quasi-geostrophic turbulence 
will be isomorphic to two-dimensional turbulence, is expected to hold. This was confirmed 
in a numerical simulation by Vallgren & Lindborg (20101, where the dissipation operator 
and the forcing term are both independent of the vertical coordinate. 

Recently, |Tulloch fc Smith ( 2009[ ) proposed a more sophisticated two-layer two- mode 
quasi-geostrophic model that has succeeded in reproducing the Nastrom-Gage spectrum. 
More importantly, using their quasi-geostrophic model, Tulloch & Smith ( 2009 1 confirmed 
that it is possible for a downscale potential enstrophy cascade to coexist with a downscale 
energy cascade, thereby vindicating the fundamental premise of the Tung-Orlando theory. 
A further advantage of the Tulloch-Smith QG model is that it seems to have resolved 
the small-scale energy dissipation sink problem implied by the HVK argument of |Smith| 
(2004). Since the Tulloch-Smith QG model is still a severely restricted approximation of 
the full quasi-geostrophic model, it is reasonable to expect that the small-scale energy 
dissipation sink problem will remain resolved under the full quasi-geostrophic model. 
This raises many interesting questions that are, unfortunately, beyond the scope of this 
paper, but will be investigated in future work. 

The goal of the present paper is to add one more piece to the overall puzzle by looking 
at the forcing range instead of the dissipation range. We will thus consider the effect of 
symmetric versus asymmetric forcing on the dynamics of the two-layer quasi-geostrophic 
model. We will first show that when the model is forced exclusively through the temper- 
ature equation, this results in antisymmetric forcing on the potential vorticity equations 
for both layers. Consequently, the energy forcing spectrum -Fe(A;) and the potential en- 
strophy forcing spectrum Fa{k) are related as Fc(k) = (fc 2 + k^)FE(k), with fc# the 
Rossby wavenumber. For forcing-range wavenumbers fc ~ fc/ <C ku, we have approxi- 
mately F G (k) ~ k 2 R F E (k). It follows that if we neglect Ekman dissipation, then the ratio 
of the enstrophy flux rj over the energy flux e injection to the downscale cascades will 
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satisfy (77/e) ~ k 2 ^, and therefore the transition wavenumber k t has to be in the vicinity 
of the Rossby wavenumber fcfj. As it turns out, this is indeed the approximate location 
of the transition wavenumber kt in the actual Nastrom-Gage spectrum as well as in the 
|Tung fc Orlando (2003a) simulation. 

We will show that asymmetric Ekman dissipation tends to decrease the ratio k t ~ 
(rj/e) 1 ^ 2 as long as the potential energy spectrum dominates the kinetic energy spec- 
trum in the forcing range. This peculiar behaviour results from the asymmetry of the 
effective forcing between the two layers, caused by the introduction of the Ekman term 
into the bottom layer. This claim is further supported by our consideration of the ran- 
dom Gaussian forcing model, in which the bottom-layer forcing is directly suppressed 
by a controlled scalar factor. Unfortunately, there are a number of open questions and 
outstanding issues concerning the distribution of energy between potential energy and 
kinetic energy. Furthermore, the modeling assumption that the Ekman term suppresses 
forcing only at the lower-layer is equivalent to ignoring the unknown effect of the entan- 
glement of the Ekman term with the interlayer interaction, and that is the underlying 
problem. 

It should be noted that, in the context of the two-layer model, unless the dissipation 
terms at small scales can dissipate the energy and potential enstrophy at the same rate 
with which they are injected to the downscale range, the downscale energy and potential 
enstrophy cascades will simply fail to develop. It is not yet obvious, in terms of theory, 
whether the two-layer quasi-geostrophic model can dissipate this much energy, a problem 
previously discussed by Gkioulekas & Tung (2007a I. On the other hand, in the real 



atmosphere, we note that at larger wavenumbers, the dynamics transitions from quasi- 



geostrophic to stratified three-dimensional turbulence. According to Lindborg (2007), 



the transition to stratified turbulence occurs at a scale of about 100km. Since stratified 
turbulence, like three-dimensional Navier-Stokes turbulence, does have an anomalous 
energy dissipation sink, it follows that any amount of energy injected at large scales can 
and will be dissipated. Furthermore, since potential enstrophy continues to be conserved 
under stratified dynamics, the two cascades can continue to coexist for scales less than 
100km. On the other hand, the two-layer model is indeed realistic at the small synoptic- 
scale wavenumbers, where the forcing takes place, so using it to explain the rates of 
energy and potential enstrophy injection at the forcing range is a fair argument. 

As we have explained above, in the Tung-Orlando theory, the location of the transition 
wavenumber k t is directly determined by the relative magnitude of the downscale poten- 
tial enstrophy flux 77 over the downscale energy flux e. A different mechanism underlies an 



SQG model that was recently proposed by Tulloch & Smith (2006). In their model, there 



is only one cascade, whose scaling exponent changes with wavenumber k, because the 
self-similar scaling of the model's nonlinear term changes with increasing k. As a result, 
in the Tulloch-Smith SQG model the transition wavenumber k t is strictly constrained to 
coincide with kn, because the transition in the scaling of the nonlinear term of the model 
occurs at kn. This dynamics of the Tulloch-Smith SQG model is analogous to that of the 



LANS a- model of three-dimensional Navier-Stokes turbulence (Holm, Livescu, Jeffery, 
Kurien, Taylor fc Wingate|2005 ), in which, once again there is a single downscale energy 



cascade with k~ b/6 scaling in the energy spectrum, with a transition to a steeper k~ 3 
slope at higher wavenumbers k, because of the introduced distortion of the Navier-Stokes 
nonlinearity. In both cases, we are essentially dealing with a single homogeneous solution, 
associated with a single flux coefficient, which, to first approximation, we can consider 
bifractal. 



Under the quasi-geostrophic models used by Tung & Orlando (2003 a I and Tulloch & 
Smith (2009), on the other hand, we are dealing with two independent homogeneous 
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solutions each of which, to first approximation, can be considered monofractal. Because 
the two solutions are independent of each other, as far as the nonlinearity is concerned, 
it is possible for the transition wavenumber fc t , in principle, to have any arbitrary value, 
since its location is determined solely by the relative magnitude of the two homogeneous 
solutions against each other. From the viewpoint of this paper, the constraint k t <~ kji 
is a weak indirect constraint that originates from the combined effect of anti-symmetric 
forcing and the large-scale Ekman dissipation term on the energy and potential enstrophy 
injection rates. This paper argues that the two-layer quasi-geostrophic model is consistent 
with placing k t near Ur. 

The paper is organized as follows. In section 2 we review the two-layer quasi-geostrophic 
model and introduce various preliminaries. In section 3 we derive the potential enstrophy 
and energy forcing spectra for the case of a generalized multi-layer quasi-geostrophic 
model, and discuss the random Gaussian forcing model. These results are applied to the 
two-layer quasi-geostrophic model itself in section 4. Conclusions and discussion are given 
in section 5 and technical matters are discussed in the appendices. 



2. The two-layer model equations 

The two-layer model is defined by two vorticity-streamfunction equations and a tem- 
perature equation which read 

^+J(^ 1 ,Ci+/) = -^w + di ) (2.1) 

d ^ + J^2X2 + .f) = + 2 [oj + d 2 + e 2 , (2.2) 

We see that the temperature T is advected by the average streamfunction (ipi + V>2)/2. 
Here -01 and 4>2 are the streamfunctions of the top and bottom layers; w is the vertical 
velocity; Ci = V 2 i/>i and Q 2 = V 2 ?/^ are the relative vorticities, and d\, d 2 , and e 2 are 
the dissipation terms given by 

di = (-1) k+ VV 2k Ci = (-ir + W 2K+ Vi, (2.4) 
d 2 = (-1) K + VV 2K C 2 = (-1) k+ VV 2k+2 V> 2 , (2.5) 

&2 = -VeC,2 = -VE^ 2 tp2- (2.6) 

The terms d\ and d 2 represent momentum dissipation of relative vorticity and e 2 rep- 
resents Ekman damping from the lower boundary layer. Furthermore, h is the height 
between the top and bottom rigid horizontal boundaries (the two vorticity layers and the 
temperature midlayer divide the space between the horizontal boundaries into four equal 
intervals, with the temperature midlayer situated between the two vorticity layers), / is 
the Coriolis term, TV is the Brunt- Valsala, frequency, and Q the thermal forcing term. 
The temperature T is related with the streamfunctions ipi and ip 2 via the geostrophic 
condition T = (2/h){ipi — ip 2 ). Finally, J(a,b) is defined as the Jacobian between the 
fields a and b and it reads: 

5) ® a ^ ^ ® a (2 7) 

' dx dy dx dy 

Solving for the vertical velocity to in the temperature equation, after substituting the 
geostrophic condition, leads to the definition of the potential vorticities qi and q 2 given 
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by 



gi = vVi + / + ^(ik-^i), (2.8) 

q 2 = ^ 2 + f-^(f 2 -^i), (2.9) 
and their corresponding governing equations which read: 

d S± + J^ uqi ) =f 1 + d 1 , (2.10) 

d ^ + J{^2,q2)=h + d 2 + e 2 . (2.11) 

Here, k R is the Rossby deformation wavenumber defined as k R = 2\/2 / /(hN), fi = 
— (l/4)fc|j/i(5o, and f 2 = {1/ A)k 2 R hQt). The derivation is shown in Appendix [Aj Although 
the argument is well-known folklore, we want to note mainly that: (a) the dissipation 
terms have the same form in the relative vorticity equations as they do in the potential 
vorticity equations; (b) the thermal forcing term Q appears on both top and bottom 
potential vorticity equations with opposite signs. Consequently, both layers are forced 
anti-symmctrically by the same forcing term, except with opposite signs. 

It is also well-known that the two-layer model, in the absence of forcing and dissipation, 
conserves the total energy E given by 



E{t) = - [^i(x J *)? 1 (x J *)+^(x,i)gi!(x > i)]dx, (2.12) 
Jr 2 

and the potential enstrophies G\ and G 2 for each layer given by: 

Gx{t)= f <z 2 (x,i)dx, (2.13) 
Jr 2 



G 2 {t) 



f <z 2 2 (x,t)dx. (2.14) 
Jr 2 



To properly define all the relevant spectra associated with these conserved quantities, 
consider first the Fourier expansions of the streamfunctions fields ip a (x, t) and potential 
vorticity fields q a (x,t) (a — 1,2): 

^ a (x,i) = / ?/> Q (k,t)exp(ik-x) dk, (2.15) 
Jr 2 

q a (x,t) = / g a (k, t) exp(ik • x) dk. (2-16) 

Jr 2 

In Fourier space, the potential vorticity fields q a and streamfunction fields ip a are related 

by 

q a (k,t) = ^L afi (\\k\\)i, a (k,t). (2.17) 



Here, the sum runs over all layers, in this case /3 = 1,2, and L a p(k) is a wavenumber 
matrix defined as 



-k 2 fcfj/2 +k 2 R /2 
+^fl/2 — k — k R /2 



(2.18) 



In real space, the same relation between the potential vorticity q a and the streamfunction 
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ip a can be written in terms of a corresponding differential operator J§f Q ^, as follows: 



t)=y)jSfa/J^«(x,t). (2.19) 







It is easy to see that the matrix L a p(k) is non-singular, for k > 0, and can therefore be 
inverted. The inverse matrix L~l(k) defines a corresponding inverse integrodifferential 



operator Jz?~p . Note that in Eq. ( 2. 17 1 we have neglected the (3 contribution to the Coriolis 
term /, since, for the case of our planet, the impact of the /^-effect on the Nastrom-Gage 
energy spectrum is negligible. We have also neglected the latitude dependence of /, on the 
premise that we are interested in the ensemble average of the energy spectrum restricted 
on a thin strip of the Earth's surface that is oriented parallel to the equator. These 
approximations cause the Coriolis term / to drop out of the nonlinear Jacobian terms 
altogether. 

Let us now introduce the following notation. Consider any arbitrary abstract scalar 
fields a(x) and 6(x), which can be snapshots in time of either the streamfunction fields 
tp a (x,t) or the potential vorticity fields q a (x,t) for a given level a. Let a <fc (x) and 
6 <fc (x) be the fields obtained from a(x) and b(x) by setting to zero, in Fourier space, 
the components corresponding to wavenumbers whose norm is greater than k. Formally, 
a <fc (x) is defined as 

a <fc (x)=/ dx / dk g(fc 7 jj ko|l) cxp(ik ■ (x-x ))a(x ) (2.20) 
Jr2 Jr2 4tH 

dx P(fc|x-x )a(x ), (2.21) 



with H (x) the Heaviside function, defined as the integral of a delta function: 

r x ( 1, if x e (0, +oo) 

H(x)= / S(t)(1t=1 1/2, ifx = . (2.22) 

Jo { 0, if x€ (-oo,0) 

Here P(fc|x — xo) is the corresponding low-pass filter kernel. Obviously, 6 <fe (x) is defined 
similarly. We now use the two filtered fields a <k (x) and b <k (x) to define the bracket 
(a,b) k as: 



(a,b) l 



d 

dk 



f dx (a <fc (x)6 <fc (x)) (2.23) 
Jr 2 

1 / dfl(A) ([a* (kAe)b(kAe) + a(kAe)b* (kAe)}) . (2.24) 

2 Ja£SO(2) \ ' 



Here, d(k) and 6(k) are the Fourier transforms of a(x) and fr(x), SO(2) is the set of all 
non-reflecting rotation matrices in two dimensions, dtt(A) is the measure of a spherical 
integral, e is a two-dimensional unit vector, and (•) represents taking an ensemble average. 



The star superscript represents taking the complex conjugate. Also note that Eq. (2.23) 



is the definition of the bracket, and Eq. (2.24) is a consequence of the definition 



It is easy to see that the bracket is both symmetric and bilinear, in that it satisfies: 

(a,b) k = (b,a) k , (2.25) 
(a,b + c) k = (a.b) k + (a,c) k , (2.26) 
(a + b,c) k = (a,c) k + (b,c) k , (2.27) 

Moreover, every (a/3)-component of the operator jSf a/ g is self-adjoint with respect to the 
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bracket, which gives 

(3? a0 a,b) k = (a,3? a0 b) k =L a0 (k)(a,b) k , (2.28) 
and the same property is also satisfied by every component of the inverse operator ■ 

(■&«pa,b) k = ( a ^«P b ) k = ■ ( 2 - 29 ) 

Using the bracket, we now define the energy spectrum E{k) = — (?/>i, qi) k — (tp2,Q2)k 
and the potential enstrophy spectra G%(k) — (?i,<ft} fe and G2(k) = {q2i<h)k for each 
layer. We also define G(k) = G\(k) + G2(k) as the total potential enstrophy spectrum. 
This method of defining spectra was previously used by Frisch ( |1995 ), and it is equivalent 
to the standard definition of spectra in terms of narrow spherical shells in Fourier space 
(see Eq. (2.24)). It is also superior in that one can generalize the definition of spectra to 
non-homogeneous flows by removing the spatial integral in Eq. (2.23), thereby obtaining 
location-dependent spectra. 

It is useful to be able to rewrite the above spectra in terms of a streamfunction spectrum 
Gap(k) = (V>ct) i>p) k - From the bilinear property of the bracket, it follows that the energy 
spectrum E(k) reads: 



= -^L a0 (k)C a0 {k). 

Likewise, the potential enstrophy spectrum G(k) reads: 

G{k) = ^ (la,Qa) k = ^2( ^2-&a/3'lpj3,^2-5? a y'<P 1 



(2.30) 
(2.31) 

(2.32) 



a/3 



^ L aj3 (k)L ai {k)C l 3 1 {k). 



(2.33) 
(2.34) 



Stated in this way, these expressions easily generalize to multiple-layer quasi-gcostrophic 
models simply by using a different matrix L a0 (k) with more rows and columns. 



3. Forcing Spectrum in general 

Now let us turn our attention to the forcing spectra of the two-layer model. We begin 
by writing the potential vorticity equations in the following more abstract form: 

7^ + J(j>a, Qa) = ^ Pg^P + fa- (3-1) 

Here, T> a p is a matrix operator that accounts for all the dissipation terms and f a is 
the forcing term acting on the a-layer. Using this abstract formulation will shorten the 
calculations below considerably. For the case of thermal forcing, the forcing terms take 
the form f\=ip and fa — —(p. 

Multiplying both sides with the inverse operator Jt?~a gives the following governing 
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equation for the streamfvmctions 

dt 



E ^JJ(^, %) = E KeVfhti + E ^U- 

13 MP 



(3.2) 



Now, let us define a streamfunction-forcing spectrum (frapik) = (fai' l Pp}k an d recall our 
definition of the streamfunction spectrum C a p(k) = (ip a ,ipp) k . Differentiating C a p{k) 
with respect to time gives 



dC a p{k) _ / dip a 



dt 



dt 



and we may write a governing equation for C a p(k) in the form 
dC aP (k) 



dt 



+ T a p(k) = D aP (k) + F a p(k). 



(3.3) 



(3.4) 



Here, T a p(k) is the contribution from the nonlinear Jacobian term, D a p(k) is the contri- 
bution from the dissipation term, and F a p(k) is the contribution from the forcing term. 
Our interest here is in the forcing contribution F a p(k). It is easy to see that F a /3(k) can 
be written in terms of the streamfunction-forcing spectrum <fi a p{k) as follows: 



F a p(k) = I E # ) + ( V-a, E *frf' 

\ 7 / fe \ 7 

= E L-^k) (/ 7 , + E L£(k) Wot 

7 7 

= Et^w^w + L-^{k)^ ia {k)]. 



k 
7/fe 



(3.5) 
(3.6) 
(3.7) 



Here, we have replaced the time derivatives in Eq. (3.3) with the forcing term from 
Eq. (3.2 1. We now note that governing equations for the energy spectrum E{k) and the 
potential enstrophy spectrum G(k) can be obtained by applying the operators indicated 
by equations (2.31 ) and (2.34) to the governing equation for the streamfunction spectrum 
C a p(k). These equations are analogous to Eq. (3.4) and they take the form: 

dE{k) 



dt 
dG{k) 
dt 



T E (k)=D E (k) + F E (k) 1 
+ T G (k) = D G {k)+F G (k). 



(3.8) 
(3.9) 



Here, the terms above have analogous definitions. The next step is to write the forcing 
spectrum F E (k) for the energy and F G (k) for the potential enstrophy in terms of F a g(k). 
Using the operator indicated by Eq. ( |2.31 ), the energy forcing spectrum F E (k) reads: 

F E (k) = -Y, L ^i.k)F aP {k) = -E^a/ 3 (fc)K 7 1 (fc)^(fc) + L^{k)4> la {k)} (3.10) 

a/3 a/37 



E 



Ewfc)^(fc) 



(*0 



I3y L \ a 



E 

ocy 



^2L^(k)L^(k))^ a (k) 



01 



(3.11) 
(3.12) 
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A similar calculation gives the potential enstrophy forcing spectrum Fc(k). We use the 
operator indicated by Eq. (2.34), and we find that Fc(k) reads: 



F G (k) = Y J L a p{k)L ai {k)F Pl {k) 



^ L a ~f(k) 



J2L a0 (k)L^(k) 



(k) + Y,L a p(k) 



a/35 



Y,L ai (k)L-£(k) 



a~/5 a(3S 



= L a p(k)4> a p(k). 



(3.13) 
(3.14) 

(t>sp{k) 

(3.15) 
(3.16) 



Note that the FE(k) calculation is dependent on the symmetry assumption L a p(k) = 
Lp a (k), which multi-layer quasi-geostrophic models do satisfy. On the contrary, the Fc{k) 
calculation is not dependent on this symmetry assumption. 

Because of the dependence of the streamfunction-forcing spectrum 4> a /3(k) on the 
streamfunction ip a , it is not possible to predict the forcing spectra Fsik) and Fc(k) 
solely from the statistical properties of the forcing term f a . The sole exception is the 
case where f a is a random Gaussian field that is delta-correlated in time such that 



(/ a (xi,ti)/ /3 (x 2 ,t 2 )) = 2Q ce/ 9(xi,x 2 )<5(*i - t 2 ), 



(3.17) 



where Q Q ^(xi, x 2 ) is assumed to be known. Then, it can be shown that the streamfunction- 
forcing spectrum (j>a/3(k) is given by 



b a p(k) = Y,Qai(k)L^(k). 



(3.18) 



Here Q Q7 (fc) is the correlation spectrum of the forcing term f a given by 



Qq,^ (fc) 



_d_ 
dfc 



J dxdydz P(fc|x - y)P(k\x - z)Q aj (y, z). 



We give a detailed derivation of this result in Appendix [Bl 



(3.19) 



4. Estimating the downscale injection rates 

We will now consider three different arguments for estimating the ratio ij/e of the 
potential enstrophy injection rate 77 to the energy injection rate e to the downscale inertial 
range. It should be noted here that a careful distinction needs to be made between the 
total injection rates to the flow, given by integrating Eq. (4.1) and Eq. \AA\ , versus 
the partial injection rates to the coexisting downscale cascades of potential enstrophy 
and energy. For the purpose of estimating the transition wavenumber k t ~ yjrjje, it is 
the partial downscale injection rates 77 and e that are relevant. These partial downscale 
injection rates are dependent on both the forcing term and the Ekman term, and that 
complicates the task ahead. 

In the first argument, we assume that the model is forced exclusively through the 
temperature equation, and we disregard the effect of the Ekman term to the downscale 
injection rates. In the second argument, we consider the combined effect of thermal forc- 
ing and the asymmetric Ekman dissipation term to the energy and potential enstrophy 
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partial downscale injections rates. We will show that the Ekman term tends to shift the 
transition wavenumber k t towards large scales, but this claim is predicated on the hy- 
pothesis that the potential energy spectrum dominates the kinetic energy spectrum at 
the forcing range, and the status of this hypothesis is presently uncertain. This prompts 
us to consider the third argument, where we force both potential vorticity equations with 
random delta-correlated in time Gaussian forcing where the forcing on the bottom layer 
is suppressed using a control factor fi. The underlying modeling assumption is that the 
asymmetric Ekman term suppresses forcing on the bottom layer but not at the top layer. 

By combining our results from these three arguments, we will argue that the two-layer 
quasi-geostrophic model tends to place the transition wavenumber k t near the Rossby 
deformation wavenumber k^. It should be noted that, due to interlayer interaction, it 
is not obvious whether the Ekman term actually suppresses predominantly the bottom- 
layer forcing. Furthermore, for models where the forcing is flow-dependent, there is further 
uncertainty on the effect of the Ekman term on forcing and the overall adjustment of the 
partial downscale injection rates. These caveats are discussed further below. 

We begin the argument by rewriting our general expressions for the energy forcing 
spectrum Fsik) and the potential enstrophy forcing spectrum Fa(k) in terms of the 
streamfunction- forcing spectrum </> Q| g(fc) for the special case of the two-layer model by 



substituting the corresponding matrix L a p(k) from Eq. (2.18). It is easy to see that the 
energy forcing spectrum Fs(k) reads: 

F E (k) = -2^0 QQ (fc) = -2[<M*0 + M*)]. (4- 1 ) 

a 

Likewise, the potential enstrophy forcing spectrum Fo(k) reads: 

F a (k) = 2'£ f L a p(k)<t> a p(k) (4.2) 

a/3 

= 2[L 11 (k)4> n {k) + L 12 {k)4> 12 {k) + L 21 [k)<t>2i{k) + L 22 (k)4> 22 {k)] (4.3) 
= -2(fc 2 + 4/2)[0 n (fc) + <j> 22 {k)} + 2(4/2)[0 12 (fc) + ^ (fc)]. (4.4) 

Using these expressions as our point of departure we now proceed with our analysis of 
the three energy and potential enstrophy partial downscale injection rate estimates. 

4.1. Estimate 1: Thermal forcing neglecting Ekman term 

Under thermal forcing we may assume that the potential vorticity equations are forced 
with /i = (p and f 2 = —(p. Let us define $i(/c) = {(p,ijji) k and $ 2 (fc) = (ip,ip 2 ) k as 
the streamfunction correlators with ip. We may therefore write the components of the 
streamfunction- forcing spectrum 4> a p(k) as: 

Mk) = Mi) h = *i{k), (4.5) 

0ia(*) = <¥>,lk> fc = *2(AO, ( 4 -6) 

M*) = <-V,^>* = -*i(*), (4-7) 

& 2 (fc) = = -* 2 (fc). (4.8) 
It follows that the energy forcing spectrum Fg{k) is given by: 

F E (k) = -2[0 n (fc) + 0a 2 (*)] = -2^) - $ 2 (fc)], (4.9) 
and the potential enstrophy forcing spectrum Fc(k) is given by: 

F G (k) = -2(fc 2 + *&/2)[&i(fc) + 22 (fc)] + 2(4/2) [(Mfc) + hi(k)} (4.10) 
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= -2(k 2 + 4/2)[$!(fc) - $ 2 (fc)] + 2(4/2)[$ 2 (fc) - $!(fc)] (4.11) 
= -2(fc 2 + 4)[$ 1 (fc)-$ 2 (A : )]. (4.12) 

We see that the energy forcing spectrum ^(fc) and the potential enstrophy forcing 
spectrum Fc(k) are related as 

F G (k) = (k 2 + k 2 R )F E (k). (4.13) 

In the limit k ~ fc/ -C k R we find that FQ(k) ~ feni^ffc), which implies that the ratio 77/e 
of injected potential enstrophy rate 77 to injected energy rate e is approximately equal 
to the square of the Rossby deformation wavenumber fc#. It follows that, if all of the 
injected energy and enstrophy cascade downscale and get successfully dissipated at small 
scales, we will then have a double potential enstrophy-energy cascade with transition 
wavenumber k t ~ k R . We suggest therefore, with some caveats to be discussed further 
below, that the two-layer quasi-geostrophic model provides a competent explanation for 
why the Nastrom-Gage spectrum exhibits a transition from k~ 3 scaling to fc~ 5 / 3 scaling 
near the Rossby deformation wavenumber k R . 

It should be stressed that this calculation neglects the effect of Ekman dissipation of the 
energy and potential enstrophy injection rates, and is only good as a first approximation. 
In the next two arguments we will suggest that the Ekman term may tend to decrease k t 
further and that it is unlikely that it can suppress the partial downscale energy injection 
rate, as one typically expects in two-dimensional turbulence. 



4.2. Estimate 2: Thermal forcing with asymmetric Ekman dissipation 

Let us now consider the effect of asymmetric Ekman dissipation on the partial downscale 
rates of potential enstrophy and energy injection. It is well-known that in two-dimensional 
turbulence, large-scale dissipation predominantly dissipates most of the injected energy 
while allowing a considerable fraction of enstrophy to cascade to small scales. This is 
evidenced by all of the numerical simulations that have successfully reproduced the en- 
strophy cascade with k~ 3 spectrum ( Kaneda k, Ishihira|2001 Lindborg fc Alvelius|2000 



Pasquero & Falkovich 2002). If the same thing were to occur in the two-layer quasi- 



geostrophic model, it would undermine our previous argument concerning the location 
of the transition wavenumber k t , because in that argument we assumed that most of the 
injected energy cascades downscale. 

As far as the downscale cascades are concerned, they "feel" forcing from both the 
forcing term and the Ekman term. It is therefore necessary to define the effective forcing 
spectra ^(fc) and S^cik) in which the effects of asymmetric forcing and Ekman dissi- 
pation are included together. We use calligraphic notation to distinguish them from the 
forcing spectra F^ik) and Fu(k) defined via Eq. (4.1 1 and Eq. (4.4). The partial injection 
rates ij and e to the downscale cascades are given by integrating ^(fe) and $c{k). As a 
matter of mathematical expediency, we can still calculate ^(fc) and ^cik) via Eq. (4.1 ) 
and Eq. (4.4) by redefining the forcing terms /1 and / 2 to include the asymmetric Ek- 
man term. The calculation below shows that the effect of the asymmetric Ekman term 
is to tend to decrease the effective energy forcing spectrum ^(fc), thereby acting as an 
energy sink. However, as long as the potential energy spectrum dominates the kinetic 
energy spectrum at the forcing range, the effective potential enstrophy forcing spectrum 
5g (fc) will decrease much faster, thereby shifting the transition wavenumber to smaller 
scales. It should be stressed that in most quasi-geostrophic models, the forcing term is 
flow-dependent, so the reaction of the flow to the effect of the asymmetric Ekman term 
adds another degree of uncertainty to the above claims. On the other hand, for random 
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Gaussian forcing that is delta-correlated in time, the reaction of the flow does not affect 
the effective forcing spectra. 

Let us now elaborate on the above argument. We begin by redefining fi = (p and 
,/ 2 = — ip — v E \7 2 ip 2 , as discussed above. Recall that C 12(h) — (ipi, ip 2 ) k , and let us define 
U\(h) = (V'ijV'i)k aim U%(k) = (tp2,ip2)k- K follows that the streamfunction-forcing 
spectra <f> a p(h) are given by: 

0ii(fc) = (V > Vi) fc = *i(*) > (4-14) 
0ia(*) = (v.^)fc = 4a(*), (4-15) 
hi(h) = {-V-v E V 2 Mi) k = ~<t>i(k) + v E k 2 C 12 (k), (4.16) 

022(*) = (-^-^V 2 V2,V>2> fe = -$ 2 (fc) + ^fc 2 L/ 2 (fc). (4.17) 

Substituting to Eq. ( |4.1[ ), we find that the effective energy forcing spectrum 3 E (k) is 
given by: 

J E (k) = -2[$i(fc) - * 2 (fc) + ^fc 2 C/ 2 (fc)] (4.18) 
= F E (k)-2u E k 2 U 2 (k). (4.19) 

Since U 2 (k) is positive-definite, we have U 2 (k) > 0, and therefore the asymmetric Ekman 
dissipation term decreases the rate of energy injection. Note that if the forcing term ip is 
depend ent on the flow, as is the case in both the Tung-O rlando and Tulloch-Smith simu- 
lations ( |Tulloch fc Smith||2009| |Tung fc Orlando||2003a[ ), the decrease by the v E k 2 U 2 (k) 
term could be counteracted by an increase by the 2[$i(fc) — $ 2 (fc)] term, as pointed out 
to the author by an anonymous referee. However, if tp is given as a standard independent 
random forcing term, which is common practice in turbulence simulations, then J E (h) 
is decreased by the Ekman term. 

Similarly, the effective potential enstrophy forcing spectrum Jc(h) is given by 

2g(*) - (fc 2 + h 2 B )3 E (k) = k 2 R [<f> n (h) + cf> 12 (h) + 21 (A;) + <fi 2 2(k)] (4.20) 

= iy E k 2 R k 2 [C 12 (k) + U 2 (k)}, (4.21) 
and it follows that 3 r G(k) is given by 

J G (fc) = (k 2 + k 2 R )3 E (k) + v E k 2 R k 2 [C 12 (k) + U 2 (k)}. (4.22) 

We see that whether the Ekman term tends to shift k t upscale or downscale depends on 
the sign of Ci 2 (A;) + U 2 (k). It is already known that U 2 (k) > 0. However, C\ 2 (k) can be 
either positive or negative. The condition C\ 2 (k) < is necessary but not sufficient in 
ensuring a transition wavenumber shift towards large scales. The necessary and sufficient 
condition for effecting such a shift is U 2 (k)+Ci 2 (k) < 0. On the other hand, the condition 
C\ 2 (k) < U 2 (k) is sufficient to ensure that the Ekman term dissipates potential enstrophy, 
which is expected on physical grounds. To show this, we rewrite the effective potential 
enstrophy forcing spectrum 3^a(k) in terms of the potential enstrophy forcing spectrum 
Fc(k) as follows: 

?a(k) = (fc 2 + k 2 R )[F E (k) - 2v E k 2 U 2 (k)\ + v E k 2 R k 2 [C 12 (k) + U 2 (k)} (4.23) 
= F G (k) + y E k 2 [k 2 R (C 12 (k) - U 2 (k)) - 2k 2 U 2 (k)}. (4.24) 

Consequently, the Ekman term dissipates potential enstrophy if k R (C\ 2 (k) — U 2 (k)) — 
2k 2 U 2 (k) < for all wavenumbers k in the forcing range. Since U 2 (k) > 0, due to being 
positive-definite, it is sufficient that C\ 2 (k) < U 2 (k). 

We can gain some insight on C\ 2 (k) by relating it with the kinetic and potential 
energy spectra E K (k) and E P (k) which are defined as follows: Let tp = (ip 1 + tp 2 )/2 and 
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t = (tpi - ip 2 )/2. So, ipi = ip + t and ip 2 = ip - r. Following [Salmon| ( 1978} |1980|), the 



definitions of the spectra E K {k), E P (k), and Ec{k) in terms of tj) and r are given by: 



E K (k) = 2k 2 (ip,ip) k , 
E P (k) = 2(k 2 + k 2 R )(T,r) k , 
E c (k)=2k 2 (V»,r) fc . 



(4.25) 
(4.26) 
(4.27) 



It can be shown that the definitions are self-consistent, i.e. E(k) — Exik) + Ep(k). It is 
easy now to write C\2(k) in terms of E^{k) and Ep{k): 



Ci 2 (k) = (ip u ip 2 ) k = (ip + T,iP-T) k = (iP,iP) k - (ip, r) k + (r, rP) k - (r, r). 



E K (k) _ E P {k) 
2k 2 2{k 2 + k 2 R ) ' 



(4.28) 
(4.29) 



We see that requiring Exik) <C Ep(k) for all wavenumbers k in the forcing range is 
sufficient to ensure that C\ 2 {k) be negative. 

To obtain a necessary and sufficient condition, we first note that 



U 2 (k) 



It follows that 



{f 2 ,i> 2 ) 

E K (k) 
2k 2 



= (ip — T,1p 

E P (k) 
2(k 2 + k 2 R ) 



E c (k) 
k 2 ■ 



{r,r 



U 2 (k) + C 12 (k) 



E K (k)-E c (k) 
k 2 



(4.30) 
(4.31) 

(4.32) 



and therefore U 2 (k) + Ci 2 (k) < if and only if Ex{k) < Ec(k) for all wavenumbers k 
in the forcing range. It should be stressed that, as far as the transition wavenumber kt is 
concerned, the relevant requirement is that the ratio Jci^/^Eik) should be decreased 
by the Ekman term. It is easy to see from Eq.(4.22) that 9 r G(fc)/3 r £;(fc) is a linear function 
of ve with slope k 2 R k 2 [C\ 2 {k) + U 2 (k)]. Thus, the condition Exik) < Ec(k) is indeed the 
needed necessary and sufficient condition. 

Without a detailed phenomenological understanding of the two-layer model, it is hard 
to say whether this condition is satisfied. Salmon (1980) has argued that in the two- 



layer model, energy is being injected as potential energy and gets converted to kinetic 
energy near the Rossby wavenumber k R . We may therefore expect the potential energy 
to remain dominant in the forcing range, provided that most kinetic energy does not 
inversely cascade back to large scales again. Nevertheless, this is an open question that 
should be carefully investigated via numerical simulations. In the next section we will 
provide an alternate argument supporting the claim of a transition wavenumber shift to 
large scales, predicated on the hypothesis that the Ekman dissipation term suppresses 
forcing only at the lower-layer, thereby assuming that the interlayer interaction does not 
propagate Ekman dissipation into the top layer. The spectrum C\ 2 {k) captures, in effect, 
an aspect of the dynamics of this interlayer interaction. 

Be that as it may, we note that it is also possible to formulate arguments that sug- 
gest the opposite conclusion as follows: As Tung & Orlando (20036) have shown, an 



equipartition of kinetic and potential energy is expected in the extreme baroclinic limit 
represented by the SQG model (i.e. Ep(k)/EK(k) — 1). The opposite limit, if generalized 



for all scales, is the case of three-dimensional stratified turbulence where Lindborg (2006) 



observed a 1/3 distribution of the total energy between potential and kinetic such that 
Ep(k)/EK(k) ~ 1/3, with the exact ratio being somewhat dependent on the rotation 
rate. For both cases we have Ex(k) > Ep(k). Furthermore, in a recent direct numerical 
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simulation of the full quasi-geostrophic model by Vallgren & Lindborg (20101, it has 



been confirmed that the total energy spectrum E{k) is equipartitioned between poten- 
tial energy Ep{k) and the two horizontal components of kinetic energy, leading to an 
approximate ratio Ep(k) / Ex{k) ~ 1/2, consistent with the theory of Charney (1971). 
A deviation seems to occur in the forcing range where Ep(k)/ Ex(k) ~ 1, but any value 
between 1/2 and 1 violates the sufficient condition Ejc(k) <C Ep(k). It should be noted, 



however, that the simulation of Vallgren & Lindborg (2010) is forced symmetrically in- 



stead of antisymmetrically, and it is uncertain how that may affect the partition ratio 
of energy between kinetic energy and potential energy. Also uncertain is the effect of 
restricting the full quasi-geostrophic model to two layers and using asymmetric instead 
of symmetric dissipation. 

At this point, one could argue that if the Ekman term in Jcik) is negligible, then it 
doesn't matter either way whether jQ{k) is increasing or decreasing. We will now argue, 
using a phenomenological order of magnitude estimate, that the Ekman adjustment of 
the potential enstrophy forcing spectrum is not expected to be negligible. The 

argument is as follows: On the assumption that most injected energy cascades to small 
scales, 3^E(k) is proportional to the downscale energy flux e. If we also assume that the 
forcing spectrum is spread over a wavenumber interval with width proportional to the 
average forcing wavenumber kf, then we get the dimensional estimate ^E(k) ~ s/kf. 
This estimate is a lower bound for 3^ E (k) since, as an anonymous referee noted, it is 
possible, in principle, for the forcing spectrum to be concentrated on a peak with width 
Ak narrower than kf. We also assume that k 2 [Ci2(k) + U2(k)], which has the dimension 
of the energy spectrum, scales as k 2 [Ci2(k) + [7 2 (fc)] ~ rj 2 / 3 kj 3 , consistent with the 
downscale potential enstrophy cascade spectrum. Putting these two phenomenological 
estimates together, for forcing-range wavenumbers k ~ kf <C kp, we estimate the two 



terms on the right-hand-side of Eq. (4.22) as: 



A = (k 2 + k 2 B )3 E (k) ~ ek 2 R /k f , 
■B = v E k 2 R k 2 {C 12 (k) + U 2 {k))~ 



(4.33) 
(4.34) 



Using the relation 77 ~ sk 2 between the potential enstrophy flux r\ and the energy flux e, 
we find that the ratio of the two terms is estimated by: 



A 



ek 2 R kf 1 



sk 2 



v k 2 



VET] 



2/3 



vet\ 



v E kpV 2/3 k f 

For the potential enstrophy flux r\ we use the value r\ 



2/3 k 2 



A/3 



v E 



(4.35) 



10" 15 s- 3 estimated by Cho fc 



Lindborg (2001) by structure function anal ysis. For the Ekman coe fficient i>e, we use 
the number ve ~ (6.7days) _1 ~ 10 _6 s _1 by Tung & Orlando (2003a). Finally, from the 
Nastrom-Gage spectrum itself, we can estimate k t ~ 10 _a km _ and kt ~ 10 _4 km _1 for 
the transition and forcing wavenumbers. Using these numbers, we find that A/"B ~ 10 _1 , 
which implies that the terms A and 23 are comparable within one order of magnitude, 
so the effect of the constant coefficients is likely to play an important role in deciding 
which term is dominant. Note that if the forcing spectrum is concentrated on a peak with 
width Afc with Ak -C kf, then ^(fc) is increased, thus the change to the ratio A/H in 
turn indicates a diminishing impact of the Ekman term on the transition wavenumber 
kt- Consequently, within the framework of the above phenomenology, our estimate of the 
A/ 15 ratio represents a worst-case scenario, in the sense that the effect of the Ekman 
term can't be stronger than this estimate. 

In light of the above, it is very important to further investigate, with numerical simula- 
tions, the effect of the Ekman term on the injection rates, using both the quasi-geostrophic 
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model of Tung & Orlando (2003 a I and the quasi-geostrophic model of Tulloch & Smith 



(2009). Specifically, for the case of the two-layer model, future numerical studies should, 



at the very least, investigate the interlayer spectrum C\2(k) and the partition of energy 
between kinetic energy and potential energy. 

4.3. Estimate 3: Asymmetric random forcing 

In the previous case, we have seen that the effect of asymmetric Ekman dissipation on 
the forcing range is to tend to decrease the rate of energy injection. There is, however, 
ambiguity regarding whether the enstrophy injection rate is increasing or decreasing, and 
whether the Ekman term shifts the transition wavenumber kt towards small scales or large 
scales. The underlying problem is that, due to the effect of the layer to layer interaction 
on the relationship between potential vorticity and streamfunction, it is not obvious 
whether dissipating the bottom layer streamfunction ip2 is equivalent to dissipating the 
bottom-layer potential vorticity 02. On the other hand, we will show now that if the 
bottom-layer forcing is directly suppressed via a control factor /i, that will indeed result 
in a reduction of the ratio Fc{k) / F E {k) in the forcing range. 

To that end, let us assume that the forcing terms for the top and bottom layers 
respectively are fi = <p and /a = — f-V with < \i < 1. Decreasing ju increases the 
suppression of the bottom-layer forcing term f 2 . We also assume that ip is a delta- 
correlated in time random Gaussian field with correlation spectrum Q(fc). In appendix 
|Cj we show that the forcing-streamfunction spectra (j) a p(k) can be expressed in terms of 
Q(fc) as follows: 

-Q(fc)lMfc) 



2fc 2 (fc 2 



Here, ipa/3 are given by: 



^n(fc) = (2k 2 + k 2 R ) - nk 2 R , 
^ 12 (k) = k R - l x(2k 2 + k 2 R ), 
^2i(fc) = -M(2A: 2 + 4)+M 2 4, 
= + ^ 2 (2fc 2 + k 2 R ). 



(4.36) 

(4.37) 
(4.38) 
(4.39) 
(4.40) 



Without explicitly calculating the forcing spectra F E (k) and i*b(fc), we can readily argue 
that since 



^n(fc) + iMfc) + V>2iW + ifai{k) = 2(k 2 + 4)(1 - M ) 2 , 



it follows that 



F G (k) - {k 2 + k 2 R )F E (k) = 4[^n(fc) + fc 2 (fc) + <j>2i(k) + <Mfc)] 

_ -Q(fc)hMfc) + iM*0 + iM*0 + fe(fc)] 



2k 2 (k 2 + k%) 
-Q(fc)[2(fc 2 + fc|)(l- M ) 2 ] 
2k 2 {k 2 + k 2 R ) 



4(1 - M ) 2 Q(fc) 



k 2 



via Eq. (4.20) and therefore 



F G (k) = (k 2 + k 2 R )F E (k) 



fc|(l~ M ) 2 Q(fc) 
k 2 



(4.41) 

(4.42) 
(4.43) 

(4.44) 
(4.45) 



We note that since the third term in the equation above is always negative, suppressing 
the lower- level forcing leads to a large-scale shift of the transition wavenumber k t . For 
fi = 1, as expected, we recover the previously derived relation Fc(k) = (k 2 + k 2 R )F E (k). 
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Another way of looking at the problem is by explicitly calculating the ratio Fq (k) / Fe (fc) 
and showing that it decreases with decreasing /i = 1. As shown in appendix [C] an explicit 
calculation of the forcing spectra Fsik) and Fq(K) gives: 

/i 2 )fc 2 + (1 - fi) 2 



F E (k) 
F G (k) 



20(fc)[2(l 

2fc 2 (fc 2 

2Q(fc)(l + /i 2 ). 



(4.46) 
(4.47) 



For the antisymmetric case fi = 1, the energy forcing spectrum Fg;(fc) reduces to 
Fs(k) = 4Q(fc) / (fc 2 + k R ) and the enstrophy forcing spectrum Fa(k) reduces to Fc{k) = 
4Q(fc), thereby recovering our previous more generally applicable result Fc(k) — (fc 2 + 
k 2 R)FE{k), which suggests a transition wavenumber k t ~ k Rl in the limit fc ~ kj <C k R . 

For the extreme case fi — 0, whereby the bottom-layer forcing is completely suppressed, 
again under the limit fc ~ fc/ -C fc_R, the energy forcing spectrum Fsik) is given by 



F E (k) 



2Q(fc)[2fc 2 



4] 



2fc 2 (fc 2 



2Q(fc)4 

01,21,2 

fv rv 



fc 2 ■ 



(4.48) 



It follows that 
-2fc/. 



and the enstrophy forcing spectrum Fc(k) is given by Fc(k) = 2Q(k). 
Fc(k) ~ 2k 2 FE(k), which suggests a reduced transition wavenumber fc t • 

The two extreme cases /i = 1 and /i = indicate that the ratio Fc(k)/ FE(k) decreases 
with decreasing /i from approximately k 2 R to 2fc 2 . We can confirm that this is indeed 
the case by taking the partial derivative with respect to the parameter fi. The partial 
derivative reads 



d 

dfi 



F G (k) 
F E {k) 



0_ 

dfi 



2(l + M 2 )Q(fc)2fc 2 (fc 2 
2Q(fc)[2(l + /i 2 )fc 2 + (l 



-4) 
rfk 2 R ] 



2k 2 {k 2 + k 2 R ) 



d_ 



2fc 2 (fc 2 + fc 2 ? ) 



[2(1 + /i 2 )fc 2 + 
2fc 2 j(l- At )(l + M ) 



[2(l + ^ 2 )fc 2 + (l 



M) 2 fc|] 2 



For < < 1, it is easy to see that every factor is positive, and therefore 



d_ 
dfi 



F G (k) 



F E (k) 



> 0. 



(4.49) 



(4.50) 



(4.51) 



(4.52) 



Consequently, the ratio Fc(k)/FE(k) decreases with decreasing fi. We conclude that if 
the asymmetric Ekman damping term on the bottom-layer streamfunction %p2 indeed 
suppresses the effective forcing of the bottom-layer potential vorticity, then the ratio 
Fc{k) / F E {k) will tend to decrease, thereby indicating a tendency to reduce the transition 
wavenumber fc t . 



5. Conclusions and Discussion 

In the present paper, we have sought out to explain why the transition from fc~ 3 
scaling to fc~ 5 / 3 scaling in the Nastrom-Gage spectrum occurs near the Rossby defor- 
mation wavenumber k R , where the atmospheric turbulence is still governed under quasi- 
geostrophic dynamics instead of three-dimensional dynamics. According to the Tung- 



Orlando theory ( Tung fc Orlando||2003« I , the entire Nastrom-Gage spectrum represents 
a downscale potential enstrophy cascade that co-exists with a downscale energy cascade. 
The location of the transition wavenumber fc t is thereby controlled by the ratio r)/e of 
the downscale potential enstrophy flux r\ over the downscale energy flux e and given 



two-layer quasi- geostrophic turbulence 



19 



by k t ~ \/r]/e. That ratio, in turn, depends on the large-scale forcing and the effect of 
large-scale dissipation on the injection of potential enstrophy and energy. 

We have shown that in the two-layer quasi-geostrophic model, which is reasonably 
applicable in the forcing scales, thermal forcing leads to antisymmetric forcing of the 
potential vorticity layer equations. This, in turn, yields a ratio r\je of the potential 
enstrophy injection rate r\ over the energy injection rate e that is approximately equal 
to k\. So, if most of the injected potential enstrophy and energy cascades towards small 
scales, then the transition wavenumber kt will be approximately equal to fcfj. 

At this point, one might object by arguing, drawing from an analogy with two- 
dimensional turbulence, that the large-scale Ekman dissipation will get rid of most of the 
injected energy at the forcing range while allowing a considerable amount of potential 
enstrophy to cascade to small scales. As it turns out, it is far from obvious that the 
two-layer model behaves in this manner. In general, the Ekman term always dissipates 
some amount of energy, and may or may not dissipate potential enstrophy, depending 
on the sign and magnitude of the interlayer spectrum C\2{k). We have shown that if the 
potential energy spectrum Ep(k) dominates the kinetic energy spectrum Ex(k) in the 
forcing range, then the downscale potential enstrophy injection rate r\ will be dampened 
faster than the downscale energy injection rate s. The resulting reduction in the r]/s 
ratio will tend to shift the transition wavenumber k t towards large scales. This tendency 
becomes exact, if the forcing used in these simulations is made independent of the flow. 

Unfortunately, there is some ambiguity in the results of our direct analysis of the Ek- 
man term, due to the dependence of the direction of the transition wavenumber shift 
on the spectral distribution of the energy between kinetic and potential energy. Using 
a random Gaussian forcing model, we have shown that, under the assumption that the 
Ekman term suppresses forcing predominantly at the bottom layer, the ratio rj/e will 
be decreased, thereby shifting the transition wavenumber kt to larger scales. While this 
assumption may seem obvious, on physical grounds, it requires us to disregard the possi- 
bility of Ekman dissipation being propagated to the top layer via the nonlinear interlayer 
interaction. Without a more detailed understanding of the phenomenology of the two- 
layer model, and especially the interlayer spectrum C\2{k), this is as far as we can go on 
this problem in terms of theory. 

Another problem with our argument is that it is only one-half of the whole story. In 
order for the injected potential enstrophy and energy to form a steady-state cascade, it 
is also necessary that the small-scale dissipation terms be able to dissipate the poten- 
tial enstrophy and energy at the same rate with which they are injected. In a strictly 
two-dimensional model, this is impossible, because the potential enstrophy and energy 
fluxes n<3(fc) and IIe(A;) are constrained by the inequality k 2 UE(k) — IIc;(fc) < 0, for 
all wavenumbers k not in the forcing range (Gkioulekas &: Tung 2005«|& ). However, as 
we have shown previously in Gkioulekas & Tung (200761, the asymmetric Ekman dis- 
sipation term can potentially cause this flux inequality to be violated. If that occurs, 
then a transition from k~ 3 to fc~ 5 / 3 scaling is possible near the wavenumber kt where 
the aforementioned flux inequality breaks down. Unfortunately, it is not easy to derive a 
rigorous necessary and sufficient condition for violating the flux inequality, in the form of 
a lower bound for v El without introducing phenomenological assumptions. In light of the 
controversy with the Tung-Orlando simulation ( |Gkioulekas fc Tung||2007a Smith|[2~004 



Tung|2004[ ), this energy dissipation sink problem remains an open question. On the other 



hand, we are quite certain that this flux inequality was successfully violated in the more 
sophisticated two-mode two-layer quasi-geostrophic model of Tulloch & Smith (2009), 
which produced coexisting cascades of potential enstrophy and energy consistent with 



20 



Eleftherios Gkioulekas 



the Tung-Orlando theory. We do not yet have a detailed mathematical understanding of 
how this violation came about. 

Ultimately, the question of whether QG models can break the flux inequality is some- 
what academic, albeit interesting. As Lindborg (2007) has shown, at scales less than 
100km, the assumptions that underlie the quasi-geostrophic model break down. This 
breakdown acts in our favor by giving us an anomalous energy dissipation sink at large 
k, thereby further facilitating the breakdown of the flux inequality. What is less obvious 
is whether there is still an effective potential enstrophy dissipation sink at small scales, 
occurring either at length scales where the flow is still stratified or via a violation of 
potential enstrophy conservation at even smaller scales where the flow becomes entirely 
three-dimensional. If yes, then we have a full accounting of the entire process: quasi- 
geostrophic dynamics is thus responsible for injecting potential enstrophy and energy 
at a proportion leading to k t ~ and three-dimensional dynamics is responsible for 
dissipating both at small scales. If no, then the widely accepted interpretation of the k~ 3 
part of the Nastrom-Gage spectrum as a downscale potential enstrophy cascade is itself 
in jeopardy, regardless of whether or not one agrees with all other aspects of the Tung- 
Orlando theory. An alternate explanation of the Nastrom-Gage spectrum as a downscale 
helicity cascade (with fc -7 / 3 scaling instead of k~ 3 ) coexisting with a downscale energy 
cascade is the only remaining hypothesis on the table, if we were to completely rule out 
quasi-geostrophic dynamics for all length scal es (|Bershadskii, Kit fc Tsinober|1993 



nover, Eidelman, Golbraikh fc Moiseev||1999 
Moiseev fc Chkhetiani||1996| ). 



Chkhctiani, Eidelman & Golbraikh 



Bra- 



2006 



It is fair to say that this paper does not resolve all of the outstanding controversies 
with respect to the Nastrom-Gage spectrum. For example, we have not yet completely 
resolved the energy dissipation sink issue in the Tung-Orlando simulation, or the question 
of whether the k~ 3 part of the Nastrom-Gage spectrum is a helicity cascade or a potential 



enstrophy cascade. In spite of extensive numerical evidence, e.g. by Koshyk, Hamilton fc 



Mahlman| fl999[) |Koshyk fc Hamilton| fl2001|), |Skamarock| ([2004D , [Takahas hi, Hamilton 



fc Ohfuchij ( |2006[ ), and Hamilton, Takahashi fc Ohfuchi| ( |2008| ) , I believe that both ques- 



tions are still open at the present time. Furthermore, within the framework of the theory 
presented in this paper, we have posed the new open question of the effect of Ekman 
dissipation on shifting the transition wavenumber away from the Rossby wavenumber 
Underlying all this, is the theoretical question of whether the location of the transition 
wavenumber k t is flexible and controlled via the magnitude of the two fluxes associated 
with two independent coexisting cascades, as proposed by Tung & Orlando (2003a), or 



whether it is inflexible and pinned down near the Rossby wavenumber kn by a scaling 
transition inherent in the nonlinearity, as typified by the Tulloch-Smith SQG model (|Tul 



loch fc Smith|2006| . While we are advocating for the flexible placement of the transition 



wavenumber k t , it is fair to say that the question deserves further scrutiny. 

An anonymous referee has also raised the question of whether gravity waves can play 



a role in the Nastrom-Gage spectrum, as was conjectured by Dewan (1979) and VanZadt 



( 1982 ). It is well-known that gravity waves vanish in the quasi-geostrophic limit, therefore 
they are not expected to be relevant over the quasi-geostrophic range of length scales, as 
were rigorously determined by Lindborg (2007). According to Gage & Nastrom (1986), 
the agreement between the measured wavenumber spectra and frequency spectra, sug- 
gests that the spectrum arises from strong turbulence and not from gravity waves. [Gage 
& Nastrom (1985) also noted that "the energy levels and shapes of the horizontal and 



vertical energy spectra are not consistent with existing models of internal wave spectra" , 
with the caveat that the inconsistency could be originating from shortcomings of these 
internal wave spectral models. Given these arguments against the gravity wave interpre- 
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tation of the Nastrom-Gage spectrum, and the "folklore" belief that quasi-geostrophic 
dynamics does not allow a downscale energy cascade, it was necessary for |Tung fc Or-| 
lando (2003a I to demonstrate that the entire Nastrom-Gage spectrum can be reproduced 



entirely by quasi-geostrophic dynamics in order to bolster their hypothesis of coexisting 
cascades of potential enstrophy and energy, even under very restricted two-dimensional 
approximations of quasi-geostrophic dynamics. Tung & Orlando (2003a) however did 



acknowledge that gravity waves could play a role in enabling the coexistence of the two 
downscale cascades. 

As for the gravity wave interpretation, many relevant questions are still not settled. 



For instance, Terasaki, Tanaka & Zagar (2011) point towards a very interesting possibil 



ity: Using reprocessed observational data provided by the Japan Meteorological Agency, 
they showed that after decomposing the total energy into a quasi-geostrophic component 
and a gravity waves component, the quasi-geostrophic component yields a fc -3 potential 
enstrophy cascade contribution spanning the entire range of resolved length scales, and 
the gravity wave component yields a A; -5 / 3 energy cascade contribution coexisting over 
the same range of scales. The total energy spectrum is thus the linear superposition of 
the two contributions. 

This picture is consistent with the Tung-Orlando theory and the linear superposition 
hypothesis p roposed by|Gkioulekas k. Tung ( 2005a|6 ) and Gkioulekas & Tung (2006). As 
explained by |Gkioulekas k, Tun g (2006), the underlying principles involved are universal 
and originate from the linearity of the underlying statistical mechanics, so we expect them 
to remain valid, beyond two-dimensional turbulence, in all related dynamical systems 
that allow the coexistence of cascades of energy and enstrophy. Under the scenario of 
coexisting quasi-geostrophic and gravity wave dynamics, indicated by |Terasaki et al.\ 
(2011), the transition wavenumber is still entirely controlled by the injection rate ratio 



77/e, given the confirmed validity of the linear superposition principle. The remaining 
open question is whether the main results of this paper concerning the injection rates 
ratio (i.e. r\je ~ fc^) can be generalized even beyond quasi-geostrophic models. We believe 
that further research is needed in that direction. 



It is a pleasure to thank Ka-Kit Tung and Joe Tribbia for discussion and correspon- 
dence. The idea of an energy-enstrophy flux inequality was originally communicated to 
me in e-mail correspondence with Sergey Danilov. All anonymous referees also provided 
very valuable feedback that went a long way into improving the paper. 



Appendix A. The potential vorticity streamfunction equations 



In this appendix, we derive the potential vorticity equations Eq. (2.10) and Eq. (2.11 ) 



from the relative vorticity equations Eq. (2.1) and Eq. (2.2) and the mid-layer tempera- 



ture equation Eq. (2.3). Our goal is to demonstrate that the potential vorticity equations 



are forced anti-symmetrically, a key property for the argument of the present paper, and 
that the dissipation terms in the relative vorticity equations retain the same form in the 
potential vorticity equations. The derivation is dependent on the following properties of 
the Jacobian J(a, b): 



J(a, b + c) — J(a. b) + J(a, c), 
J(a + b,c) = J(a, c) + J(b, c), 
J(a,a) = and J(a,b) = —J(b,a), 
J(Aa,/i6) = A/zJ(a, b), 



(Al) 
(A 2) 
(A3) 
(A 4) 
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where A and \x are constants. 



The first step is to solve for the vertical velocity u> in the temperature equation Eq. ( 2.3 ) 



From the geostrophic constraint T — (2//i)(0i — ^2) we write the advection term in the 
temperature equation as: 

]{ i>i^h T \ = i[jty 1) fa)-jfa t fo) + J fo ) fa)-jfo i1 fo)] (a 5) 

= -fj(^,0 2 ). (A6) 



It follows that the vertical velocity ui reads 

/ 

2/ 



Qo 



and therefore 



hN 2 



h 2 



wOi - ih) - J(ipw>p2) - 



J^(01 - 02) - J(lpl,lp2) ~ 



(A7) 
(A 8) 

(A 9) 



Here we have defined the Rossby deformation wavenumber fc# = 2\/2//(WV). 

The next step is to define the potential vorticities q± and qi for the top and bottom 
layers correspondingly as: 



?i = vVi + / + ^(ik-M, 

92 = V 2 2 +/~ -^(^2- 0l)- 



(A 10) 
(All) 



The advection terms J(ipi,q±) and J(ip2,q2) of the potential vorticities with respect to 
the streamfunctions 0i and 02 are given by: 



and 



J(V'i.gi) = J(AXi + f + (4/2)(V>2 _ ^)) 

= J(^,Cl+/) + ^J(01,0 2 ), 



J(0 2 , «a) = J(*h, Ca + / - (4/2)(02 - 0i)) 



(A 12) 
(A 13) 



(A 14) 
(A 15) 

(A 16) 

Therefore, differentiating the top-layer potential vorticity q\ with respect to the time t 
gives: 

ggi 
dt 



dt 



k% d 



(02 



2 dt 

-J(sl>i,b + f) 
-J(0i,Ci + /) + 



-01) 

2/ 



w + di + 



2 di v ^ 2 w 



d 



— (01 - -02) - J(01,0 2 ) - 



k\d_ 
2 dt 



(A 17) 
(A 18) 

(02 - 0l) 

(A 19) 
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J(ipi,qi) -Q + dx, 



k 2 hk 2 
= -J(^,Ci + /) - -fJ{^2) - ~^Qo + di 
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with Q defined as Q — hk^Qo/A. Likewise, differentiating the bottom-layer potential 
vorticity q2 with respect to the time t gives: 



dq 2 
dt 



d(2 
dt 



fdt^ 



-JW>2,C2 + /)- 
-^2,C2 + /)- 

-^2,C2 + /)- 
-J(ip2,q 2 )+Q- 



2/ 



ft 
2 



e 2 



^2 

2 <9r^ 2 



V2) - J(lpl,1p2) 



hQ 



- d 2 + e 2 . 



e 2 



(A 22) 
(A 23) 
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(A 24) 

(A 25) 
(A 26) 



The governing equations Eq. (2.101 and Eq. (2.111 for the potential vorticity follow. 



Appendix B. Streamfunction-forcing spectrum under random 
Gaussian forcing 

Let us consider the case of a generalized multi-layer model forced at each layer a with 
random Gaussian forcing f a such that 



(fa(x-i,ti)fp{x2,h)) = 2Q a/ 3(x 1 ,x 2 )S(t 1 - t 2 ). 



(Bl) 



From the Novikov-Furutsu theorem (Furutsu 1963 Novikov 1965) we know that, given 



a functional R[f], the correlation between f a and R[f] reads 



(f a (xi,ti)R\f\) 



dx 9 



dt 2 (/ Q (xi,ti)/ (x2,t 2 )) 



8R[f] 



5fp(x 2 ,t 2 ) 



(B2) 



It should be noted that implied is a space-time approach in which it is the entire forc- 
ing history / that is being mapped to a number by the functional R[f]. The ensemble 
average is understood to average over all possible forcing histories. The idea is to treat 
the streamfunction tp a of layer a at a given point in space-time as a functional of the 
entire forcing history, and then use the Novikov-Furutsu theorem to evaluate the forcing- 
streamfunction spectrum 4> a p(k). This idea follows a similar argument by |McComb| ( |1990| ) 
for the three-dimensional Navier-Stokes equation. The argument proceeds as follows: 
Recall first the definition of the filtering kernel: 



a <fc (x,i) = / P(fc|x-y)a(y,t)dy. 
By definition, the streamfunction-forcing spectrum <f) a p(k) is given by 

4> a p(k) = (f a ,i>p) k = ( /«>5Z ) = H L i3lW (/«' ?7)fc 

\ 7 1 k 1 



(B3) 

(B4) 
(B5) 
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dxdydz P(fc|x - y)P(fc|x - z) (/« (y, i)g 7 (z, *)) 



2\3 



(Bl 



Using the Novikov-Furutsu theorem, we calculate the forcing-streamfvmction correlation 
(/a(y, i)g 7 (z, t)}, and find that it reads: 



ifa(y,t)q 1 (z,t)) 



<%yQM) 
<5q 7 (z,t) 



dw / dt (fa(y,t)f s (w,t )} 

[ dw f dt 2Q aS (y,w)5{t-t ) 
Jr 2 Jr \oM w >*oJ 

2/ dwQ tt5 (y,w)/^ (M) 

JR 2 



<*/i(w,t) 



(B7) 
(B8) 
(B9) 



To evaluate the variational derivative of potential vorticity g 7 (z, t) with respect to layer 
forcing fg(w,to) we first note that, by causality, the potential vorticity q 1 (z,t) at time 
t is related with the initial potential vorticity g 7 (z, 0) at time to = by an equation of 
the form 

q 7 (z,t) = g 7 (z,0) + / d*o>T 7 [g(to)](a)+ / dt / 7 (z,t ). (BIO) 
Jo Jo 

Here, N 7 [g(io)](z) represents the combined effect of the nonlinear and dissipation terms. 

The third integral represents the causal contribution of the forcing term. Let us assume 

now that < r < t, and differentiate the above equation variationally with respect to 

fg(w,T). We immediately find that 



5q 1 {z 1 t) 



$fs(w,T) Sfg(vf,T) 

= A r j(z, t; w, t) + 23 7<5 (z, t; w, r 

with A^g(z, t; w, r) and 23 7< 5(z, t; w, r) given by 

(5 



dt X r [«(* )](a)+ / dt Uz,t ) 



A-ys(z,t;w,T) 



<*/«(w,t) Jo 
5 

5f s {-w,T 
S 



dt ^[q(t )]{z) 



dt Q f 7 (z,t ) 



dt 



S^[q(t )](z) 
*/«(w,r) 



/ dto / dz 5(z - z )/ 7 (z ,to) 
Jo JR 2 

/ dt / dz - t )(5(z - z )/ 7 (z ,to) 

JR JR 2 



5 7 aiZ"(r - i)<J(z - w). 



(Bll) 
(B12) 

(B13) 
(B14) 
(B15) 

(B16) 
(B17) 



Here, H(t) is the previously defined Heaviside function. For Eq. (B13), we rely on the 
principle of causality to restrict the integral from r to t. It is easy to see that for t = r, 
Ay<s(z, i; w, r) and ® 7 a(z,t;w, r) simplify to: 



■A 7 a(z,i; w,i) = J dt 



£N 7 [g(*o)](' 



= 0, 



S 7 5 (z, t; w, t) = -5 7 5^(z - w), 



(B18) 
(B 19) 
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and therefore the variational derivative of g 7 (z,i) with respect to /a(w, t) is given by 



<5q 7 (z,i) 1 
Sf s (w,t) = 2 



<5 7(5 <5(z-w). (B20) 



Substituting this result to Eq. ( B 9 ) , we show that the forcing-streamfunction correlation 
is given by 



(fa(y, t)q y (z, t)) = 2 / dw Q aS (y, w) ( 

JR 2 \ 



(B21) 



<y/«(w,t) 

2f dwQ a{ (y,w)Jy(z-w) (B22) 

/ dw Q tt7 (y, w)(5(z - w) = Q a7 (y,z). (B23) 

Jr 2 



Consequently, the forcing-streamfunction spectrum a ^(fc) reads 

d 
'dfc 



a/J (A) = $>^(*)^r /// dxdydz P(fc|x-y)P(%-z)(/ a (y,i)g 7 (z,t)) 



(B 24) 

= E L M^a%j jj m3 dxdydz P(fc|x ~ y)P(A;|x _ z )^(y ; z ) ( B 2 5) 

= ^Q Q7 (A;)^ 7 1 (A ; ). (B26) 

7 

The integral above defines the forcing correlation spectrum Q Q7 (fc), given by: 

Q Q7 (fc) = ^ J fj^ 3 dxdydz P(fc|x - y)P(fc|x - z)Q Q7 (y, z). (B 27) 

Our final result for the forcing-streamfunction spectrum is: 

<Mfc) = ^Q Q7 (fc)^ 7 1 (fc). (B28) 



Appendix C. The random Gaussian forcing model 

Let us consider the case of the two-layer quasigeostrophic model forced with f\ = if at 
the top layer and / 2 = —[Atp at the bottom layer. Here, /i is a suppression constant with 
< ji < 1 and ip is a random Gaussian field that is delta-correlated in time such that: 

(<p(xi,h)<p(x2,t 2 )) = 2Q(x 1 ,x 2 )<5(t 1 -i 2 ). (CI) 

From Q(x 1 ,x 2 ) we define the corresponding correlation spectrum Q(k) as: 

Q(k) = ^J fj^ 3 dxdydz P(fc|x - y)P(fc|x - z)Q(y, z). (C 2) 

It follows that, for a, ft € {1,2}, f a and fp are correlated according to 

(/a(xi,*l)/fl(X2,*2)>=2Q afl (x 1 ,X 2 )5(t 1 -t 2 ), (C3) 

with the components of Q a p given by 

Qu(xx,x 2 ) =Q(x 1 ,x 2 ), (C4) 
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Qi 2 (xi,x 2 ) = Q 2 i(xi,x 2 ) = -/xQ(xi,x 2 ), 
Q 22 (xi,x 2 ) = ^ 2 Q(xi,x 2 ). 

The spectrum Q a p(k) of Q Q ^(xi,x 2 ) is defined as 
d 



dk 



2)3 



dxdydz P(k\x - y)P(fc|x - z)Q a p(y, z), 



consequently its components read: 

Qu(fc) = Q(fe) ) 

Qi 2 (k) = Q 21 (k) = - f iQ(k), 

Q 22 (k) = fi 2 Q(k). 



(C5) 
(C6) 



(C7) 



(C8) 
(C9) 
(CIO) 



In Appendix [5] we have shown that under general delta-correlated in time random Gaus- 
sian forcing, the general form of the forcing-streamfunction spectrum 4>ap{k) reads: 



b aP {k) = Y,^{k)L-^{k). 



(Cll) 



We would now like to reduce this result to the case of the two-layer quasi-geostrophic 
model. Starting from Eq. (2.18), a simple calculation shows that the inverse matrix 

L ap( k ) is s iven b y ; 



L~Uk) = 



2k 2 (k 2 + k 2 R ) 



k 2 



k 2 
2/c ~\~ k 



(C12) 



We note that the inverse matrix L~l(k) is defined for all wavenumbers k > 0. Combin- 
ing the above two equations we find that the components of the streamfunction-forcing 
spectrum 4> a p{k) are: 



n (fc) = Qu{k)L^{k) + Q 12 (k)L u \k) = Q(fc)[L n 1 (fc) - ^{k)\ 
_ -Q(k)[2k 2 + k%- iik%] 



2k 2 (k 2 + k%) 

012 (k) = Qu(k)L^(k) + Q 12 (fc)L 22 1 (fc) = Q(k)[L^(k) - »L^(k)] 
_ -Q(k)[k% -fi{2k 2 + k 2 R )] 



2k 2 (k 2 + k 2 R ) 

021 {k) = Q 2 i(k)L n 1 (k) + Q 22 (k)L u \k) = Q(fc)[- M L n 1 (^) + M^uW] 

= -Q(k)[-^(2k 2 + k 2 R ) + ^k 2 R ] 
2k 2 {k 2 + k 2 R ) 

022 (k) = Q 21 (k)L 2 i(k) + Q 22 (k)L^(k) = Q{k)[-^{k) + » 2 L 2 i (k)\ 

= -o(fc)[-^4 + M 2 (2fc 2 + 4)] 

2k 2 (k 2 +k 2 R ) 

We may therefore write the streamfunction-forcing spectra as: 

-Q(*)lM*0 



<Paf3{k) 



2k 2 (k 2 + k 2 R ) : 



with ip a /3 given by 



^ii(A) = (2k 2 + k 2 R ) - fxk 2 R , 



(C13) 
(C14) 
(C15) 
(C16) 
(C17) 
(C18) 
(C19) 
(C20) 

(C21) 
(C22) 
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^ 12 (fc) = 4^ M ( 2 fc 2 + 4), 

ip2i(k) = -fi(2k 2 + k 2 R )+n 2 k 2 R , 

■4>22(k) = -/J,k% + fi 2 (2k 2 + k%). 
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(C23) 
(C24) 
(C25) 



From the streamfunction- forcing spectra (f)ap(k) we calculate both the potential enstro- 
phy forcing spectrum Fc(k) and the energy forcing spectrum F E (k) using Eq. (4.1 1 and 
Eq. (4.4 1. An easy calculation gives: 



(C26) 
(C27) 



^n(fc) + ^22 (fc) = (2k 2 + k R ) - nk% - (xk R + [i 2 (2k 2 + k 2 R ) 
= 2(l + ^)fc 2 + 4(l- / i) 2 , 

and therefore, for F E (k) we find that 

F E (k) = -2[0 n (fc) + </> 22 (fc)] = 2 fc 2 (fc 2 + fc|) 

= 2Q(fc)[2(l + M 2 )fc 2 + (l- M ) 2 fc 2 ? ] 
2k 2 (k 2 + k 2 R ) 

For the potential enstrophy forcing spectrum Fc(k), we use a slightly more subtle argu- 
ment, and we have: 



(C28) 
(C29) 



F G (k) = 2j2L a p(k)<j> aP (k) = 2j2L aP (k) 



a/3 



2£ 



J2L-^(k)Q aj (k) 



J2L aP (k)L^(k) 



cx-y 



= 2^Q aa (fc) = 2[Q 11 (fc) + Q 22 (fc)] 

OL 

= 2(l + {i 2 )Q(k). 



(C30) 

(C31) 

(C32) 
(C33) 



It is worth noting that the potential enstrophy forcing spectrum Fo(k) is independent 
of the matrix L a p(k) as long as L a p(k) is non-singular. The energy forcing spectrum 
F E (k), on the other hand, is dependent on the inverse matrix L~l(k). Eq. (C29) and 
Eq. ( C 33 ) are the main results of this appendix. 
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